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Abstract 

We use the Euler-Maclaurin formula and the Feynman-Kac formula to ex- 
tend our previous method of computation of the spectral action based on 
the Poisson summation formula. We show how to compute directly the spec- 
tral action for the general case of Robertson-Walker metrics. We check the 
terms of the expansion up to ag against the known universal formulas of 
Gilkey and compute the expansion up to aio using our direct method. 
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1. Introduction 

In a previous paper ^13j we devised a method based on the Poisson sum- 
mation formula to compute the spectral action up to arbitrary order for 
geometric spaces possessing a lot of symmetries such as products of spheres 
by tori. It is natural to wonder whether this method extends to compute 
the spectral action for the general case of Robertson-Walker metrics. We 
show in this paper that this is the case. We start by showing in fj2] that for 
even dimensional spheres one can use the Euler-Maclaurin formula instead 
of the Poisson summation formula to compute all terms of the asymptotic 
expansion. Moreover this method is fairly general and applies as soon as the 
spectral multiplicity is a polynomial and one has a polynomial expression 
for the squared eigenvalues. To cover the general case of Robertson- Walker 
metrics we first reduce in f|3] the computation of the spectrum of the Dirac 
operator to a one dimensional problem and then use the Feynman-Kac for- 
mula to compute the spectral action. We first prepare the ground in §^and 
[5l by performing the computation of the terms up to ag using the universal 
formulas of Gilkey. As another preparation we use, in ^the semiclassical 
approximation to get the accurate value of the first term gq. In f|7]we explain 
how to use the Feynman-Kac formula to compute all terms of the expan- 
sion, up to surface terms, using the Euler-Maclaurin formula. We check 
the even terms up to og against the universal formulas. This check is quite 
revealing because of the high degree of complexity of the universal formula 
of Gilkey for the ag term. But it leads us to a puzzle since the method 
yields non vanishing odd terms such as ai, 03 etc. in the expansion. In ^ 
we resolve this puzzle by showing that the operator used in our computa- 
tion was only similar to the Dirac operator and the natural symmetry of 
the latter entails the vanishing of the local formulas for the odd terms while 
justifying the computation of the even ones. We also show that this natural 
symmetry allows one to use the Poisson summation formula instead of the 
Euler-Maclaurin formula, and this simplification allows us to compute the 
local formula up to oio while giving an efficient algorithm to compute terms 
of arbitrary order. Finally we give in ^some open questions concerning the 
general terms of the expansion and the conceptual meaning of the use of the 
Brownian bridge. This paper is a mathematical preparation for the goal of 
exploring the cosmological implications of the spectral action principle (|16j. 

m, m, m)- 

2. S^, Euler-Maclaurin formula and 

In this section we explain the relevance of the Euler Maclaurin formula in 
the computation of the spectral action for spaces whose Dirac spectrum has 
two properties 

• The eigenvalues of the Dirac operator D are of the form zizh{k), 
where h{x)'^ is a polynomial function of x. 

• The multiplicity of ib/i(/c) is a polynomial P{k) in k. 



In [13] we considered the case of odd spheres with their round metric and in 
this case the polynomial P is even so that one can use the Poisson summation 
formula to compute the spectral action. We also treated the case of the even 
sphere S"^ using the zeta function which is straightforward to compute, and 
given by 

(1) TT{\Dr) = ^{c{s-3)-as-i)) . 

This function has a value at s = given using 

C -3 - C(-i = + — = • 

^ ^ 120 12 120 

The same computation gives the full asymptotic expansion of the trace of 

the heat kernel 9{t) = Trace(e~*^ ) for the sphere 5^, since 



Cs4s) = / Trace(e-*^') = — — / 9{t) f'^-^ 



dt 



so that (^g4{—2k) = {—l)^k\bk where the 6^ are the coefficients of the expan- 
sion 

m 

9{t) = Trace(e-*^') = Ur'^ - r^) + V but'' + 0(^+1) 

fc=0 

Now recall that the value of the Riemann zeta function at odd negative 
integers is given by 

in terms of the Bernoulli numbers. Thus one gets for the sphere the heat 
kernel coefficients 



(-1)^ , (B2k+2 



2fc+4 



(2) hk = —C^i{-2k) 

The Bernoulli numbers are the value at « = of the Bernoulli polynomials 
Bn{u) which are defined inductively as follows 



Bo{x) = l, BM = nBn-i{x), [ B^ 

Jo 



{x)dx = 0. 



Equivalently, these polynomials can be introduced using the generating func- 
tion 

oo 

" (U) 



(3) nu,t) = -r-^ = Y.B. 



n=0 



n! 



The presence of the Bernoulli numbers in ([2]) can be traced back to the Euler 
Maclaurin formula 

j^gik) = fg{x)dx + a{a)+ 9{b) ^ ^(,0-1) (5) _ gU-i)ia)) - 



k=a " " j=2 



i?m = M- / g^'^\x)Bm{x-[x])dx 



m\ 

which we shah now use to compute the expansion in the even sphere case. 
For 5"^, the eigenvalues of Dirac are the itfc with multiphcity — k). This 
shows that the eigenvalues and ±1 are not there (multiplicity zero) and 
the multiplicity of the eigenvalue k"^ of is |(A;'^ — k). We first take the 
test function 

f{u) = e-'^' 

We apply the Euler Maclaurin formula for a = 0, b = oo and 

g(x) = ^{x^ - x)f{x) = ^{x^ - x)e"*^' 



,2 



The Taylor expansion of g( at is given by 

, oo A A 

(4) g{x) = - Y,{x' - x){-irex'yn\ = --x + -x^{l +t) + ... 



3 
and one has 

g^'\o) = -^, ^^ = ^(i + t), g^'-'Ho) = Oit) , Vj>4 

Thus the only contribution of the term 

i=2 ^' j=2 

to the constant term in t is of the form 

^ B., ^ fB2 Bi\ A f I _ I \ 4 11 

3 120 



^ ^ ' 3 V 2 4 y 3 Vi2 i2oy 

J— 2 



The term ^(^)^^(^) is zero since 5^(0) = 5^(00) = 0. The integral 

b A poo 9 

g{x)dx = - {x^- x)e-'^'dx = -{t'^ - t'^) 
3 Jo 3 

does not contribute to the constant term in t. 

Let us now use ^ and the Euler Maclaurin formula to compute the full 
asymptotic expansion in powers of t. One has 



(2m -1)! ' ' 3V("^-2)! (m-1)! 
The relevant contribution is 

75 R /I / +m—2 j.m—1 

-^g^'-'Ko) = -(-i)-|^| ( + 

2m] 2m 6 \[m — 2)\ [m — l)\ 

Thus we get that the coefficient of t^'/kl is given by 



(-1) 



3\2k + 2 2k + 4j 
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The Euler-Maclaurin formula holds for any smooth function and the only 
thing one needs to care about is to control the remainder which is of the 
form 

Rm = M- / 9'-"'Hx)Bm{x - [X])dx 



m 







and whose size is controlled by an estimate of the form 

/•oo 

(5) <2C(m)(2vr)— / 



Jo 

In the case of the spectral action for we apply the Euler-Maclaurin 
formula to the function g{x) = f{tx'^){x^ — x) and the derivatives ^^"^^(0) of 
this odd function at x = are of the form, 

9'iO) = -/(O) 

5(3) (0) = 6/(0)-6t/'(0) 

5(5) (0) = 120t/'(0) -60t2/"(0) 

g(^) (0) = 2520 t^f'iO) - 840 t^f^^'^ (0) 

g(9)(0) = 60480 t3/(3)(o) -15120 tV^^^(O) 

5(ii)(0) = 1663200 tV^^^O) -332640 tV^^^O) 

5(i3)(o) = 51891840 tV^^Ho) - 8648640 tV^'^^(O) 
The Euler-Maclaurin formula then applies directly and gives the expansion 

(6) ^Tvacei f{tD^)) = T fitx^){x^ - x)dx + ^^^^^^ ^^^'^^^^ 



4 7o 120 2520 

4ir(0)t^ 31/(3) (0)^3 10331/(4) (0)^4 3421/(5) (0)t5 ^ 
10080 15840 8648640 3931200 +■■■+ 

where the remainder Rm is under control. To see this one can use ([5]) and 
the fact that 

<;(™)(^) = t(-3)/2p^(n,t,/('=)(n2)) 

where Pmiu, t, f^''\u^)) is a polynomial in u,t and finitely many derivatives 
of / evaluated at u^. In particular this gives Rm = 0{t^"^~'^^^'^) when / is 
smooth with compact support or more generally belongs to the Schwartz 
space. This implies that when / is a cutoff function, i. e. a smooth function 
with compact support with all higher derivatives f^^\0) = for all A; > 0, 
the full asymptotic expansion is of the form (up to an overall | ) 

^TVace(/(tZ)2)) = ^ f{u)ndu - 1 f{u)dn + + ) . 

The advantages of this method are that not only one controls the size of the 
remainder for arbitrary test functions / but also that it can be used generally 
for spaces for which one has a polynomial expression Q{k) for the squared 



eigenvalues of the Dirac operator and a polynomial multiplicity P{k). One 
applies the same method as above with 

g{x) = f{tQ{x))P{x). 

3. Dirac operator for Robertson- Walker metrics 

The starting point in the evaluation of a spectral action is the Dirac oper- 
ator of the relevant geometry. The (Euclidean) Robertson- Walker metric of 
dimension 4 with the symmetry of the round sphere is read from the line 

element 

(7) ds^ = df + (t) (dx^ + sin^ x {dO'^ + sin^ edip"^)) . 

The computation is most efficiently done by using a frame 6^ so that the 
line element simplifies to ds'^ = 6"'9°' where the flat indices a,b, - ■ ■ , are 
contracted with the Euclidean metric Sab- Thus the frame is given by 

(8) 0° = dt 

(9) e^ = ait)dx 

(10) 9^ = a (t) sin xd9 

(11) 6"^ = a(i)sinxsin6'd(^ 
which satisfies 



(12) 









(13) 


de^ = 


a 




(14) 


de'^ = 


a 


a 


(15) 


de^ = 


-^°A( 

a 


a 



asmx 

where a' denotes time derivative. The spin-connection is defined by 

(16) = u"^ A 9'' 

(17) w"''= ra;^"'' 
The non-vanishing components of ojabc are 



(18) OJioi = W202 = ^4^303 = — 



a 



cotx cot I 

i-Lyj ^^212 — '^313 — ) <^323 



a a sin x 

This shows that the Dirac operator 

(20) ^ = «5^ + i7'^w"'' 



/ N i\ I d 3a' \ 1 ^ 
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where the gamma matrices 7" are antihermitian satisfying (7")^ = — 1, and 

d_ 

sin X sin 9 dip 

is directly related to the Dirac operator on which is given by 



(22) £'3 = 7' i^ + cotx) +7'-^ + +7^^^ 
\ox J smx \o(f 2 J smxs 



3 



■ if 9 A . 2 1 M 1 .^ . 

Dirac53 =ia — + cot x + { ^ + n^^otO ] +ia . . . 

\ox J smx \oO 2 ) sin % sm i) dip 

where the cr^ are the Pauli matrices. In fact more precisely one has 

(23) 7°L»3 ~ Dirac^a -Dirac^s 

using the representation 

Moreover 7°-D3 has the same square as D3 since 7*^ anti commutes with L>3 
and has square —1. Thus the square of the Dirac operator is 

One can view the vectors in the Hilbert space of spinors as functions ^ (t) of 
t with values in 

L2(53,S)eL2(5-\S) 
where E is the spinor bundle on and the square norm of ^ is given by 



/ 



corresponding to the volume form for the Robertson- Walker metric 

■y/g = a'^ sin^ X sin 9 

The spectrum of the Dirac operator for the round sphere S"^ of unit radius 
is given by 

Spec(L>) = {±(^ + fc) I A: G Z, A: > 0} 

where the multiplicity of (^ + /c) is equal to 2'2l (^'''^~"^) . For d = 3 wc 
get the multiplicity {k + l)(fc + 2) for the eigenvalues (| + fc), which means 
the multiplicity — \ for all odd half-integers. Thus decomposing the 
vector valued function ^(t) in terms of the eigenfunctions of the operators 
±Dirac53 for the eigenvalue A, one reduces the problem to the direct sum of 
the one dimensional problems corresponding to the scalar operators in one 
dimension given by 

X ( Q 3a' \^ 1,9 a' , , , ,3 , 



One can move back to the standard measure dt in the variable t by the 
transformation 

(27) v{t) = a~^u{t) 

which gives the direct sum of the operators where 

fd' {n + l)\{n + l)a'\ 



(28) 



dt"^ a? a? 



which occurs with multipUcity 2(n+l)(n+2). We shall in ^and ^compute 
the spectral functions of the direct sum of the operators Hn = with the 
multiplicity given by 

(29) ^(n) =4(n + l)(n + 2) 

The obtained result will be shown to differ by a surface term from the 
correct result for the original four-dimensional eigenvalue problem, which 
corresponds to the direct sum of the operators © H~ with multiplicity 
^//(n). This important nuance coming from the sign in front of the term 
(n + I) a' will be dealt with at the beginning of fusing time reversal. 
As a simple example we take the case of the 4-sphere where a{t) = sin(t) 
and t G [0,7r]. We find that the spectrum of the operator is given by 

SpecH^ = {k^\k>n + 2} 

We can expand individually the terms Tr (e^''^") in the infinite sum: 

n+l -, n+1 

2 ' 2^ 2 ^ ^ 



Tr(e-^^")= e-'""" = -Y,e-'''" — Y^e-'''" ^ — Y^e 



k>n+2 



{n+l) + \{n + l){2 + 3n + n^)s + .. 



2^ ^ 2' 3^ 2' 

This expansion gives a divergent sum for the summation over n and is useless 
so it is better to use the zeta function. Taking f{x) = instead of 

f{x) = e~^^ in the above formula we get, using 

Yin + l){n + 2) = ^ 


the correct expression 

7i>0 n>0 k>n+2 k=2 



4. Spectral action for Robertson- Walker metrics 

In this section we give the spectral action for the Robertson- Walker metric 
starting from the universal formula of Gilkey. To do this we will have to 
specify all geometric quantities. Starting from 

(30) de" = uj"^ A 9'' 

(31) R'^ = duj"'' - uj"" A uj^^ 

so that the torsion 

(32) T" = de" - uj^^B^ = 

and dT" = -R'^O^. This is consistent with our definition of the Lorentz 
covariant derivative 

(33) Dij, = dfj, + ^Uf^ablab 

so that 

(34) [D^,D,] = ^R^,ab7ab 

(35) R^uab = dfj,U!i,ab — ^imc^vcb — < > ^ 

We also have the relations 

, ,ab ac, ,cab 

^ab ^-^abcdyQC ^ qcL 

~ 2 
0" = eldx" 

Covariant derivatives with respect to curved indices are related to those of 
tangent space indices through the equation 

(36) d^e^a - ^iiab^ub - {g) epa = 

which makes the covariant derivative with respect to tangent space vectors 

(37) Dfj.Va = dfj.Va- OJ^abVb 

equivalent to covariant derivatives with respect to vectors over the manifold 

(38) V^V, = d^V,-TP^,Vp 
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where Vi, = e"V^. Going back to the Robertson- Walker metric we get the 
components of the curvatures are 

a 

a 
-.11 



i?03 = A 

a 

i?12 = 7i A 6 



a? 



i?i3 = ^ — A ^3 



a' 



which give the relations 



a" 



(39) -Roioi — -R0202 — -R0303 — — 

a 

(40) -R1212 = -R1313 = -R2323 = 5 • 

The components of the Ricci tensor are then 

a" 

(41) -Roo = RocOc = 3— 

, , a" (o''-l) 

(42) Rii = Ricic = h 2 5 = R22 = R33 

and finally 

43 = = 6 — + ^ . 

The Weyl tensor defined by 

Cafccd = Rabcd—-^ {^acRbd — hcR-ad — ^adRbc + ^bdRac) + -^ {^ac^bd — ^adhc) R 

vanishes for the Robertson- Walker metric: 

(44) Coioi = i?oioi- J(i?oo + ^ii) + ^i? = 

Z D 

(45) C1212 = i?i2i2 - \ {Rii + R22) + = 

The Gauss-Bonnet term R*R* = le^^P^ea^^sRui^""^ Rpa'^^ gives 

24(-l + a-(0^)a--(0 
= ^(tp 



11 



For an operator of the form 

(46) D'' = - (g'^^'d^d, + Af'd^ + B) 
the first three Seeley-de Witt coefficients are 

(47) ao = j^ld''x^Tr{l) 

(48) "2 = ^/ d'x^gTv (e - li?) 

= 16vr^x 360 / ^'"^^^ + " ^^^'^ + 

(49) +60n£; + 180£;2 - 6QRE + 300^^) 

where 

(50) io^ = ^g,A-^'' + n 

(51) E = B-g^'^ {d,,uj,. + u^uj, - TP^.ujp) 

(52) il^j, = (9^a;,^ - + [u^, uoy] 
For the above simple Dirac operator wc have 

(53) OJ^ = l^fjab 

(54) E = -R 

4 



,2 

fjLvpa 



(55) J^;.. = ^i?^/S„6 



1 
4 

We thus get the first two coefficients 
and 

(57) «^ = T^x4x|^«'|xT5 = T(^- + S^j 

Next in order to compute the coefficient 04 we need the geometric invariants 

(58) ^,^, = 12 _ +L_^ + _L_^ 



//\ 2 



(a--l)^ 



(59) i?a6cdi?a6cd = 12 |^(^— J + ^4 J 

This gives 

a" fa'2 - 1) 

(60) 5i?2 - 8i?a6i?„6 - IRabcdRabcd = 264 ^ 



12 



While the surface term comes from —12\I\R + GOOE = SOR with 

UR = R-aa = Da^aR' = R" — ^aaoR' 

fa(^) a'a(3) {a"f a'^a" a"\ 
\ a a-^ a-^ I 

so that 3Tr(l)ni? gives a factor 3 X 4 X 6 = 72 in front of the parenthesis. 
This gives the final result for 04 : 

IGtt^ X 360 J \ \ a a-^ a-^ I 

(62) 

+^a" (a'2 -lj^=±-jdt (3a2a(4) + Qaa'a^^) + 3aa"2 - 4a'2a" - 5a") . 

We note that this is a linear combination of the total derivative Di? and the 

Gauss-Bonnet topological combination. Thus 04 does not contribute to the 
dynamics of the spectral action. The above formula fits with the general 
formula 

a4 = Y^^l d^x^ (-18C2^p, + IIR*R* + 12nii) 
since the Weyl tensor vanishes. 

5. Calculation of ag 

To calculate the next term in the spectral action we start with Gilkey's 
formula. This is fairly complicated, we thus partition the expression into 
four blocks. The first block depends quadratically on covariant derivatives 
of the curvature tensor and its contractions: 

^^-^^ j d^X^Tr {-ISD'^R + 17 R.,aR.,a - 2Rab;cRab;c - ^Rah;cRbc;a 

+^Rabcd;eRabcd;e + ^SRDR — SRab^Rab + 24RabRbc;ac + I'^Rabcd^Rabcd) 

The second block depends cubically on the curvature Rabcd and its contrac- 
tions but not its derivatives: 

77—9 ^ Z7 / d^'a^v^Tr (— 35i?^ + A2RRabRab — ^'^RRabcdRabcd + '^OSRabRbcRca 

167r^ X 9 X 7! y ^ 

-192RabRcdRachd + ^^RabRacdeRbcde " ^'^RabcdRabefRcdef " ^ORahcdRaecf Rbedf) 

The third block depends on the curvatures Qab a^id their derivatives as well 
as the curvatures Rabcd and their contractions: 

167r^ X 360 / '^^^^'^ {8^ab;Ab;c + ^^ab;b^ac;c + l2nabOnab 

— 12Q,ab^bc^ca — ^Ra^cd^ab^cd + ^Rab^ac^bc ~ ^R^ab^ab) 
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The fourth block depends on E, ^ab, Rabcd and its contractions 

X 360 / ^^^V^'^ (^°'^ + eOEDE + SQE,aE,a + QOE^^ + SQEnab^ab - lORBE 

—ARabE-ab — 12R.aE-a — SORE^ — 12nRE + bR^E — 2RabRabE + 2RabcdRabcdE) ■ 

In our case E = ^R.l and 0,^1, = \R^ivablab- Calculating the traces and 
using the identities 

Tr {-iablcd) = -4 {5achd - Sadhc) 
Tr {-Jablcdlef) = 4 [( 5bc (SaeSdf " ^afSde) - a i > b) - C < > d] 

we obtain the following table of conversion: 

^ab ~^ ~h-^abcd ' ^abOQab — ^ —\Rabcd^Rabcd , ^a6;c ~^ ~'k^^ibcd:e 

^ac^bc ^ ~2^acdeRbcde ; ^ab;b^ac;c ^ ~ ^Rabcd;bRaecd;e ) ^ab^cd ^ ~2^abefRcdef 

^ab^bc^ca ^ 2RabcdRaecf Rbedf 

Note that all other terms have an additional factor Tr(l) = 4 in the com- 
putation of their trace. This enables us to list all terms that appear in 
ae : 

1 /" / 1 1 1 1 



167r2 



1 12 2 

7^Rabcd;eRabcd;e + l^R^R ~ -^Rab^Rab + TTT^ RabRbc;ac 

252 360 315 105 



1 „ 1 1 „ 7 ,,, 

^abcd^abcd 



-Rabcd^Rabcd + „,-„^ ^'^ " -TT^RRabRab " .^^^ RRabcdRa 



140 2592 540 4320 

52 16 1 101 

"I" '^^^RabRbcRca RabRcdRacbd J^Q RabRacdeRbcde ~^ '^^^^^RabcdRabefRcdef 

(63) 

1 1 109 

~'^^RabR;ab ^gg Rabcd,bRaecd;e ~^ '^^'^^RabcdRaecf Rbedf 

We see that some of these terms are not so straightforward to evaluate. We 
note the following 

(64) R,a = SaoR' 

(65) Rab;c = GcRab — ^cadRdb - ^cbdRad 

thus since the Rab are dependent only on t, and uoab = 0, we then have 

(66) Rab;0 = R'ab 

a' 

(67) -Roi;l = — {RqO — Rll) = Ro2;2 = Ro3;3 

Notice that this last relation shows that one must be very careful in working 
with covariant derivatives which arise because of the split between space 
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and time coordinates. Complications start with second derivatives: 

= Cc {ecRab — ^cadRdb — ^cbdRad) — ^ccd {c-dRab — (^daeReb — ^dbeRae) 
— ^cad {ecRdb — ^cdeReb — ^cbeRde) — ^cbd {^cRad — ^caeRed — ^cdeRae) 

More complicated is 

^Rabcd = De {DeRabcd) 

= Ce {DeRabcd) — l^eefDfRabcd — ^eafDeRebcd 
— UebfDeRafcd — l^ecfDeRabfd — ^edfDeRabcf 

where 

DeRabcd = ^eRabcd — l^eafRfbcd — ^ebfRafcd — l^ecfRabfd — ^edfRabcf 

We note that all the components of the tensors such as 

R ) R;a ) R,ab > Rab j Rab,c > Rab,cd > Rabcd j Rabcd,e > Rabcd;ef 

are only functions of the variable t and not of the other variables ^! 
even though the variables X) ^ appear explicitly in some components of the 
connection Uabc- We omit the lengthy steps, and we just list the final ex- 
pressions. We start by the scalar curvature and the fact that the Laplacian 
on functions f{t) of t alone is given by □/ = ^° ^ait) 

(68) n^R =(r"+ 3-R') +3- ( R" + 3-R') 

\ a / OS V ^ / 

We also have easily R-aR;a = {R')'^- For the next two terms we use 

(69) Rab;cRab;c = {Kof + 3 {R'nf + 6 (^^^ {Roo " Ruf 

2 f a'\'^ a' 

(70) Rab;cRac;b = {R'oo) ~'~ ( j ("^00 " -?^ll) + Q — R'u {RqO — Rll) 

Then for the term in Rabc^.^ we get 



Rabcd;eRabcd;e — 12 



(i?Ul)' + (^'1212)' + 4 (^^) (i?0101 - i?1212) 



The term in RDR is straightforward. For the next term in Rab^Rab we use 

(71) Rab^Rab = RqO^RoO + Si^nDi^n 

where moreover 

(72) Moo = i?oo + 3 (^^ j i?oo + 6 (^^ j {Ru - Roo) 
and 

(73) □i?ii=i?';i + 3(^^^i?ii+2(^^^ {Roo-Ru) = aR22 = aR33 
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The next term is in RabRbc;ac and is given by 

RabRhc;ac =RoO ^-^00 + 3— (^^OO ~ -^ll) ~ ^ ^'o} ^^^^ ~ ^n)^ 

(74) + 3Rn ( + {Roo - Ru) + ^Ko^ 

Next for the term involving the full Riemann tensor Rabcd^Rabcd one has 

(75) Rabcd^Rabcd = 12 [i?010ini?0101 + i?1212nit:i212] 

where 

(76) DEoioi = <ioi + 3 i^oioi + 4 (^^^ (-R1212 - i?oioi) 

= 0-^0202 = CH-R^OSOS 

and 



(77) □-R1212 = -R1212 + 3 -^1212 + 4 (-f^oioi — -R1212) 

= Di^isis = □i?2323 
The term in i?'^ is straightforward. One has for the next term 

(78) RRabRab = R {Roo + 3-Rll) 

Similarly for the next term, RRabcdRabcd involving the full Riemann tensor, 
we get 

(79) RRabcdRabcd = 12-R (-RqIOI + -^1212) 

We now evaluate the three next terms (and remember that R22 = Rii) 

(80) RabRbcRca = RqO + 3-Rii 

(81) RabRcdRacbd = 6i?00^11^0101 + 6i?iii?22^1212 

(82) RabRacdeRbcde = 6 (-RqO + -f^ll) -f^OlOl + 12-^11-^1212 

The next term is cubic in the Riemann tensor and given by RabcdRabefRcdef- 
One gets 

(83) RabcdRabefRcdef == 24 (i?oi01 + -^1212) 

For the next one RabR-ab one gets 



a' 



(84) RabR;ab = RoqR" + S-RuR' 

a 

We still have one more term involving first derivatives of the Riemann tensor 

/ ^/ \ 2 

(85) Rabcd;bRaecd;e = 6 ( -RqIOI + 2— (-^0101 — -R1212) 1 

Finally the last term is 

(86) RabcdRaecfRbedf = 6 (-R1212 + 3-^1212-^0101 ) 
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We 



are now ready to collect all the terms in one expression 
8 J [420 a J a \ a J J 1008 ^ ^ 

- ^ ([Roof + 3 (^)' (^00 - i?ii)' + e^R'u {Roo - Rn?j 

^ r / , N 9 1 



1 
21 



+ 



105 



(^0101)^ + (^'1212)^ + 4 ^— ^ (-R0101 - ^1212)^ 

1 / // a' A 

+ R{ R +3-R'] 

360 V « / 

— -^00 ^-^00 + 3~-^oo + (-^11 ~ -^00)^ 

+3i?ii ( Rii + 3—R'ii + 2^ {Rqq — ) 
\ a y 

-^00 ^-Rqo + 3— {R'oo - R'u) - 3 ^— ^ (-f 

+3i2ii ( ( h 3^ ) (i2oo — Ru) H -Rq 

\ \ a / a 



{Roo — R 



3 

35 



+ 



^-Roioi ^-^0101 + 3— -Roioi + (-R1212 — -Roioi)^ 

+-R1212 ^-^1212 + 3— -R'1212 + 4-^ (-R0101 — -^1212)^ 

11 7 
■ 2592^' " l40^ ^^0° + ^^'1) ~ 1320^^^ + ^'^12) 

"2835 3-Rfi) - (-R00-R11-R0101 + -R11-R1212) 

— ((i?oo + -Rii) -Rqioi + 2-Rii-Rf2i2) + 22680 ('^oioi + -^1212) 

- ^ (^RooR" + 3— -Rii-R'^ — — ^-Rqioi + 2— (-R0101 — -^1212)^ 
109 .0 „„ „ 



jgQQ (-^1212 + 3-R1212-R0101) 

As an elementary test of this we evaluate the coefficient for a sphere 5^ 
where a = sin t, which implies 



which implies 
R = —12, -Roioi = -R1212 = — 1, 



^00 = -^11 = —3 
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which gives the coefficient to be 



31 
2520 



and thus 



31 7 ...... 31 



as = / sin tdt 



2520 J 1890 



which agrees with the expansion given in 

Substituting the expressions of R, Rqq, Ru, Rqiqi, R1212 into the above for- 
mula gives 

(87) 

^ a'{tra"{t) _ a'{t)^a"{t) a"(tf a'{tfa"{tf _ l_ 3 
240a(f)2 84a(t)2 ^ 120a(t) ^ 21a(t) 9o" ^ ^ 

+ 240a(t) + 84«(t) " 2o" ^'^^ ^'^^ " 1680 " 24o" 

This expression is quite elaborate and it will be a major test for our direct 
method of computation which we develop below, to check that we obtain 
the same result for og. This will be done in ^ The above equality together 
with ()56p . (j57p . and the computation of 04 determine the spectral action 

(88) I = TTf {D^/Ml,) = M^.Uao + M^j2a2 + foa^ + M-2/-2a6 
up to order -r-p- where 

/•OO /'OO 

(89) U= uf (n) d^z, /2 = / / (n) dn, /o = / (0) , /_2 = -/' (0) 

JO JO 

As 04 does not contribute to the dynamical equation for a{t), ag will have 
non-trivial effects on the standard model of cosmology. 

6. Semi-classical approximation and ao 

As shown in ^ the square of the Dirac operator for the Robertson- Walker 
metric associated to the function a{t) is intimately related to the direct sum 
of the operators if„ given by 

(f 

(90) Hn = -^ + Vn{t) 

(91) ^"(*) = ^^((" + ^)-"' 



where Hn occurs with the multiplicity /x(n) = 4(n-|-l)(n-|-2) for all integers 
n > 0. This means that the spectral action corresponds to the following 
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sum 



(92) TV/(D2)~^^(n)Tr/(i?„) 



n>0 



We take the functions f{u) = e and we look at the expansion when 
s ^ 0. At the semi-classical level the operator is 

and in this section we shall investigate what happens if we just use (for 
f{u) = e~'^") the first approximation to Tr/ (Hn) given by the semi-classical 
expression 

where the integration is over p G M and t. One has 



27r 7 V^2TrJ 



-sVr 



and thus the semi-classical approximation is given locally in t by the sum 



n>0 2V^n>0 



We can use the Euler-Maclaurin formula to evaluate the summation in n 
and the first contribution is an integral. As a function of n we are dealing 
with 

;u(n)e-*^"W = 4(n + l)(n + 2)e'^^(" -^"+2)) = hs{n + ^) 

where 

/i,(a;) = (4x2 - l)e-'^+'l^ 

Next one has 

/•OO 1 

(93) / (4^2 - l)e"(^-^)^dx = ^--^e4(-3+2'')« ^ 

2 

(^2(3 + ^))^/^ + e3(-3+^)'"V^(2+ (-l + 62)n) (^1 + Erf ^(-3 + 6)v^ 

where the error function is 

Erf(z) = ^ / e"" dv 
Jo 

and its expansion near is 
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Thus taking u = sa^'^, b = a' , one obtains the expansion of the coefficient 
of dt after multiphcation by which gives up to O [s] ^^"^ , 

a' a^a' a[-2 + 3{a'f) - 3a' + 4{a'f ("4 + 5 (af 



2s2 ^s3/2 8s Vl^/l^^fs 64a 

One needs to add the other terms of the Euler-Maclaurin formula, which are 

i=2 

The derivative Il's^ is of the order of where a is the largest integer a < 
a/2. Thus for the expansion up to 0[s]^/2 one can take only the terms 
corresponding to 

i'.4)-t'4"(f)^K'-l^"i"<l> 

(3) 

since h), is of order s. The computation of these terms gives, after multi- 



plication by 



\3 



2^^s 

which cancels the term independent of a' in +'^" ) ^nd thus the 

asymptotic expansion from the semi-classical formula gives up to 0[s]^/2, 



a 



3 ^2^1 a [-2 + 3 {a'Y) ^ _2,a' + A{a'f ^ ("4 + 5 (a' 



2s2 ^s3/2 85 12V^Vs 64a 

This formula is not nice at all since 

• It only involves a and a' and hence cannot reproduce the curvature 
terms. 

• It has non vanishing terms in odd powers of ^/s. 

• The coefficient of l/-v/s is not a total derivative. 

This is of course to be expected since the approximation uses the wrong 
assumption that p and the function a{t) commute. This assumption only 
holds in the case where a is constant, in that case one drops all terms in a' 
and only two terms remain: 



a^ a 



2s2 4s 

which agrees with [13j thm. 3. In general the only term we can trust is the 

3 

leading term ao = ^ which fits with the general formula 

Tr(l) f ^ ,4 

since Tr(l) = 4 and the volume of the three sphere of radius a is 271^0^. 
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7. FeYNMAN KAC formula and up to 06 WITH SURFACE TERMS 

In this section we use the Feynman Kac formula to compute the spectral 
action. We use the following formula (see |21] Theorem 6.6 page 54) for the 
local expression of the trace of e~^^", 



exp ^-2s j V{t + V2sa{u))di?jD[a] 



where a is the Brownian bridge (see [21j Definition page 40) which is a 
Gaussian random variable with covariance given by 

E{a{u)a{v)) = u{l — v) 

for u < V. In our case the potential V is given by —^df + V= \Hn 

V{t) = ^Vn{t) 

and the above formula becomes 

e-^^"(t,t) = -l=jexp (^-s Vn{t + V2^ a{v)) dv^ D[a] 
The technique then is that we evaluate the sum with multiplicity included, 

oo 

j;/i(n)e-^^"(t,t) 

n=0 

by using the Euler-Maclaurin formula, which means replacing the discrete 
index n by the continuous variable x = n + ^ and the sum over n by 

ks{x)dx + -ks\-] + — + ... 

3 ' 2 1 2 / 12 720 30240 

2 ^ ^ 

where the function kg is given by 

ks{x) = (4x2 - 1)-^ [ e^(''->D[a] 

2y/TTS J 



Jo 



+ \/2s a{v)) dv 



1-1 

ub = s I a a^"^ (t + V2s a{v)) dv 





where 



and 



so that 



To obtain these terms we have used x = n + | and 

—s / Vn{t + V2s a{v)) dv = u{b — x)x 
Jo 



aa~\t + V2sa{v))dv/ [ / a'-^it + V2s a{v)) dv 
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We use the Taylor expansion of a^^ and of a'a~'^ for the expressions of u and 
h to get an asymptotic expansion when s — )• 0. The terms of the expansion 
are coming from the general formula 



F{t + V2i a{v)) dv = Fit) + ^ ^^^^ 

where ^ 

Xk{a) = / a{v)^ dv 
Jo 

We get the following expansions 



s 2 (^/2a'(t)) xi(a)s3/2 2 {-3a' [tf + a{t)a" (t)) X2{a 



2 



and 

b = a'{t) + \/2a"(t)xi(a)Vs 
Uxi{afa'{t)a"{t) - 4x2{a)a' {t)a" (t) + a(t)x2(a)a(3) (t)) s ^ 

a{t) 

We are first dealing with the integral term in the Euler-Maclaurin formula. 
Thus we can use the formula ()93p . The leading term is unchanged and given 
by a(t)^/(2s^) as above. When we use the additional terms for u and b we 
find that the term in ^ in the expansion gets new terms depending on a, 
and becomes 

a{t)(-2 + 3a'{tf) 15a(t)a'(tf , ,0 3a(t) (-3a'(i)2 + a(i)a"(i)) , , 

8i + Ys + 2^ 

Performing the Gaussian integral one gets 

y xx{afl){a\ = ^ , y X2{(x)I){a\ = ^ 
and the term in - thus gives 

s ^ 

ajt) {-2 + 3a'{t)^) I5a{t)a'{t)^ 1 3a{t) {-3a' jt)^ + a{t)a" jt)) 1 
8s 2i 12 2^ 6 

_ a{t) {-l + a'{tf + a{t)a"{t)) 
~ Ts 
This agrees with the spectral action which gives for the ^ term 

while 

/a" o'2 1 
R = 6{ — + ^ - ^ 

\ a 

which is negative for the sphere and gives (using |5f | = 27r^a^) 
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We now look at the 03 term, including the corrections it gives 

ag = i (^a'(t) (-3 + Aa{t)a"{t)) + 2a{tfa^^\t) 

and one finds that it is the total derivative of the following expression which 
vanishes at both ends of the time interval 

+ ia(t)2a"(t) 
2 3 ^ ^ ^ ' 

The next term is the 04 term, one performs the computation in the same 
way as above and obtains 

04 = (- (5 + 4a' {tf) a"{t) + 3a{t)a"{tf + 3a{t) (3a'{t)a^^'> {t) + a{t)a^^'> [t] 

This term agrees with (j62p and is very interesting because it does not vanish 
for the sphere case where it gives 

llsin^f 4 11 

a4(sphere) = , a4(sphere)dt = 3 

but it is a "topological" term since it is the derivative of the following ex- 
pression 

-—a'{t) - ^-^^ + —a{t)a'{t)a"{t) + —aitfa'^^Xt) 

In this expression the last two terms vanish at the end points of the time 
interval, but not the first two. In fact the variation of the first two across 
the interval will not be zero in general except when a'{t) vanishes at the 
boundary which is the case of x . In the sphere case the derivative 
at the end points is ±1 and this corresponds to smoothing out the conical 
singularity at the boundary. 

The computation of is more complicated. To do the computation one 
needs to compute the integrals of polynomials in the Xj{a) under the Gauss- 
ian measure D[a] in order to obtain the coefficients. We list in the appendix 
the table of the integrals which are needed to compute up to oio- It gives 

a'{tfa"(t) a'{tfa"{t) a"{tf a'{tfa"{tf 1 , 3 a'{t)a^^\t) 
~ 240a(i)2 84a(t)2 ^120a(i)^ 21a(t) 9o" ^ 240a(t) 

a-(t)V3)(t) _^ (3) «(t)a(3)(t)2 _J_ J_ 

+ Ma{t) 20^^^^ ^' 1680 240 " 120^^ ^' 

+ 840 +i4o"WaW« W+ 

One checks that it agrees with the computation (j87p using the Gilkey uni- 
versal formula. It is important to see why the discrete terms in the Euler- 
Maclaurin formula do not contribute to a^. In fact the direct computation 
gives (up to a factor of 4) 

99xi(a)a'(t) 2'^^/2xi{a)a' {tf 29xi{a)a"{t) 
20V2a(t)3 15a(t)3 ^ Ib^aitf 
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and such terms disappear after the integration in D{a). In the general 
case the terms oe and ag, since they have the denominator a{t) to some 
power, will be singular at i = 0, but this singularity vanishes provided that 
a^^)(0) = in the case of og. These conditions are verified for the sphere. 
In fact the smoothness of the metric at the pole i = is equivalent to the 
condition 

a'(0) = ±l, a(2")(0) = 0, Vn. 

The results of this section leave us with a puzzle since we get the correct 
values for the even terms but the method yields non vanishing odd terms 
such as oi, 03 etc. They give surface terms whose integral vanishes but they 
seem to contradict the vanishing of the local trace for the odd terms. 

8. Full Dirac operator, Poisson summation and up to aio 

In this section we resolve the above puzzle by showing that the operator used 
above, namely the direct sum of the operators Hn = with multiplicity 
admits the same even terms in its spectral expansion as the Dirac 
operator and the natural symmetry of the latter entails the vanishing of 
the local formulas for the odd terms while justifying the computation of the 
even ones. We also show that this natural symmetry allows one to use the 
Poisson summation formula instead of the Euler-Maclaurin formula, and 
this simplification allows us to compute the local formula up to aio while 
giving an algorithm to compute terms of arbitrary order. 
We have seen in ^that the square of the Dirac operator is 

Now observe that the operator can be viewed as the "time reversal" of 
= 

(94) ^^ = -§2+^^^'^ 

(95) ^"^(*) = ^^((" + 0+«') 

and this implies that the spectral asymptotics for the direct sum of the H^^ 
with multiplicity n{n) are obtained from the above ones simply by replacing 
the derivatives a^'^^ by (— l)'^a^'^^. Thus no new computation is needed for the 
first terms and one checks that, for the Dirac operator, the non vanishing 
odd terms such as oi, 03 etc. cancel out, while the even terms remain 
unchanged. We shall in fact go further and show how to use this added 
symmetry to simplify the above computations trading the Euler-Maclaurin 
formula for the Poisson summation formula. We consider the function 



fs{x) = {x' - -)e"(^--)- 
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and the relevant sum is now 

oo ^ 

(96) E/«("+2) 

— oo 

Using the Poisson summation formula it is very well approximated by 

There is an overall factor of 2 coming from spinors and one multiplies by 
^ ] — and then integrates in Dial as above with 

u = s I a^^ {t + V2s a{v)) dv 



Jo 

and ^ 

ub = s a' a'"^ {t + V2s a{v)) dv 
Jo 

so that 

6 = ^ a'a''^{t + V2^a{v))dv/ (^j^ a'"^ {t + V2^ a{v)) dv 
Thus the relevant expression is 

(97) i f '-^S^±i-l±3ADla] 

One then repeats the same computation as above and finds that all the terms 
Qj for j odd vanish, while they are unchanged for even j. For instance when 
we compute 04 we find the following coefficients as a polynomial expression 
in the x{j) 

..U _. 315a' (f)-^ 
■^\^) 4a(t) 

x(l)2x(2) ^ -^Iglf + ^-fa'itfa"it) 
^ + ^ + la'{tYa"{t) + la{t)a"{tf 

x{l)x{^) ^ 55^1)1 _ 45a'(t)2a"(t) + ba{t)a' {t)a^^\t) 
^(2)2 ^ 1^ _ ^^a'{tfa"{t) + f a(t)a"(t)2 

^(2) ^ ^ - - ^ - f + !a(i)a'(t)a(=^)(t) 

^(4) ^ -^^^ + 18a'(t)2a"(t) - la{t)a"{tf - 3a(t)o'(t)a(3)(t) + ia(t)2a(^) 

16a{t) 64a(t) 

We then replace the monomials in the x{j) by their integral under D[a\ com- 
puted using integration by parts under the Gaussian measure and gathered 
in the table of the Appendix. 

The main result then is that this technique gives the full expansion of the 
spectral action to arbitrary order. The proof is the same as in [13] using the 
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fact that the remainder in the Poisson summation formula has flat Taylor 
expansion. 

For the og term one obtains the following expression 

as = TTT-Pfa) 

^ 10080a(t)4 ^ ' 

P{a) = 108a'(i)6a"(i) - lQSa{t)a' {tfa^"^) {t) + a' {t)\21a" (t) - bSSa{t)a"{tf 
+ ma{tfa^^\t)) - 3a(t)a'(t)3((9 - 256a(t)a"(t))a(3) (t) + ^a{tfa^^\t)) + 
3a(t)a'(t)2(-29a"(t)2+267a(t)a"(t)3-104a(t)2a"(t)a(4)(t)+a(t)(-72a(t)a(3)(t)2+ 
5a(^)(t)+3a(t)2a(6)(t)))+a(t)2a'(i)(-819a(t)a"(t)2a(3)(t)+6a"(t)(17a(3)(t) + 
13a(t)2a(5)(f))+a(t)(132a(i)a(3)(t)o(4)(t)-6a(5)(t)-5a(i)2a(7)(f)))+a(f)2(43a"(t)3- 
114a(t)a"(t)4 + 69a(t)2a"(t)2a(4)(t) + 2a(t)a"(t)(45a(t)a(3)(t)2 _ l8aW(t) + 
a(f)2a(6) (t))+a(t) (-24a(3) (t)2+l3a(t)2a(4) (t)2+16a(t)2a(3) (t)a(5) (t)+3a(t)a(6) (t)- 
o(t)3a(8)(t))) 

and for the sphere S*^ it simplifies to 

41 sin'^ t 4 41 

aafsphere) = , / asfsphere)^^ = - x 

^ ^ ^ 10080 Jo 3 10080 

which agrees with the direct computation dG]) of ^ 
With this method we can now reach aio which is given by 



aio = ^^^^^^{-n700a'{tfa"it) + 11700a(t)a'(t)7a(3)(t) 

+ 3a'(t)6(5a"(t)(-165 + 5096a(t)a"(t)) - 20A6a{ty a^^\t)) 

+ 3a{t)a'{tf{{825 - 34628a(t)a"(t))a(3) (t) + 746a{t)^a^^\t)) 

+ 3a(t)a'(t)4(3476a"(t)2 - 54054a(t)a"(t)3 + lU40a{tfa"{t)a(^\t) 

+ a(t)(10448a(t)a(3)(t)2 _ 429aW(f) - 2na{t)^a'-^\t))) 

+ 3a(t)2a'(t)3(78902a(t)a"(t)2a(3)(t)-2a"(t)(2222a(3)(t) + 2127a(t)2a(5)(f)) 

+ a(t)(-7992a(t)a(3)(t)a(4)(f) + 154a(5)(t) + 55a(i)2a(7)(f))) 

+ a(t)2a'(t)2(-11880a"(t)3 + 111378a(t)a"(t)4 - 68mAa{tf a" {tf a^'^'^ {t) 

+ 33a(t)(113a(3) (t)2 + UAait^a^^^ {tf + 192a(f)2a(3) {t)a^^) (t) - 4a(t)a(6) (t)) 

+3o(t)a"(t)(-31973a(t)a(3)(t)2+i738a(4)(t)+968o(t)2a(6)(t))-43a(t)4a(s)(t)) 

+ a(t)3a'(t)(-117600a(t)a"(t)3a(3)(t) + 66a"(t)2(211a(3)(t) + 172a(t)2a(5)(f)) 

- 2a(i)a"(t)(-19701a(t)a(3) (t)aW (t) + 693a(5) (t) + 238a(t)2a(^) (t)) 

+ a(t)(-2640a(3) (t)a(4) (i) - 2a(t)2 (778a(4) (t)a(5) (t) + 537a(3) (t)a(6) (t)) 

+ 33a(t)(271a(3)(t)3 + a(7)(t)) + 18a{tfa(^\t))) + a(t)3(2354a"(t)4 

-3a(t)a"(t)(3446a"(t)4+1243a(3)(t)2+924a"(t)a(4)(t))+66a(t)2(331a"(t)2a(3)(f)2 

+160a"(t)3o(4) (t)+6a(4) (i)2+9a(3) (i)a(5) (f)+4a"(i)a(6) (t))-a(t)3 (331a(3) (t)2a(4) (t) 

+ 482a"(t)o(4)(t)2 + 1110a" (t)a(3) (f)a(5) (t) + 448a" (t)2a(6)(t) + lla(8)(t)) 

-3a(t)4(74a(5)(t)2 + llQa^^) {t)a^<i) (t) + A9a^^\t)a^'^\t) + 6a" {t)a^^\t)) 

+ 3a(t)5a(io)(t))) 
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When we evaluate this expression for a{t) = sm(t) we obtain for the sphere 
S*^ the result 

, , , 31 sin3 1 r , , , , 4 31 
«io(sphere) = , as{sphere)dt = 3 x ^ 

which agrees with the direct computation ([6]) of ^ There is another very 
useful test of the above computations which is to compute the divergence at 
t = 0. The test is that this divergence should cancel under the hypothesis 
that a^^") = and a'(0) = 1. For instance for the coefficient oio above, one 
finds that the numerator admits, under the hypothesis that a'(0) = 1 and 

= for n < 4, the expansion 
i^a(io)(o)t8 + |2(2625a(3)(0)5 - 12712a(3)(0)3a(5)(0) + 6804a(3)(0)a(5)(0)2 
+ 2592a(3)(0)2a(7)(0) - 1152a(5) (0)a(7) (0) - 240a(3) (0)a(9) (O))^^ + O(t)io 

which is of order 9 in t provided a^^'^)(0) = and thus compensates the 
denominator a{t)^ in the formula for the local coefficient of the volume form 
a{tfdt. 

9. Conclusion 

We have developed in this paper a direct method for the computation of 
the spectral action for Robertson- Walker metrics in Euclidean formulation. 
We have shown that our method, based on the Euler-Maclaurin formula 
combined with the Feynman-Kac formula, gives the same result for the local 
terms of the expansion as the Gilkey formulas up to ag. We have computed 
the full expansion up to aio and checked its accuracy by performing concrete 
tests such as the direct comparison in the case of the 4-sphere. We have 
shown also how to control the remainder in the asymptotic expansion which 
only gives an approximate expression of the spectral action neglecting the 
role of the instanton contributions. In fact besides giving the above explicit 
formulas our results suggest the following questions: 

• Check the agreement between the above formulas for ag and aio and the 
universal formulas of [l], [2] and [22] . 

• Show that the term a2n of the asymptotic expansion of the spectral action 
for Robertson- Walker metrics is of the form P„(a, . . . , a(^"'')/a^"~^ where 
Pn is a polynomial with rational coefficients, and compute Pn- 

• Is there a conceptual physical meaning for the role of the integration over 
the Brownian bridge as a fluctuation of the time variable in passing from 
the semiclassical approximation to the full spectral action? 

10. Appendix 

We give in this appendix the table of integrals of polynomials in the Xj{a) 
under the Gaussian measure D[a] corresponding to the Brownian bridge. 
The computations are straightforward since they consist in integrating by 
parts under a Gaussian but time consuming so that it is useful to get the 
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assistance of a computer. We list below the integrals needed to compute up 
to aio- The meaning of an arrow such as xfxl — >■ is 

/ \2 / \2 7-.r 1 659 

where a{t) is the Brownian bridge and 

Xj;{a) = / a{v)''dv 
Jo 



1 - 

4 

X|X2 - 

2 3 
Xl-t-2 ■ 

xi ^ 

xjx3 - 
Xix|x3 - 

13 ' 

X2X3 
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